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ABSTRACT 


Aims. We study the effects of letting dark matter and gas in the Universe couple to the scalar field of the symmetron 
model, a modified gravity theory, with varying coupling strength. We also search for a way to distinguish between 
universal and non-universal couplings in observations. 

Methods. The research is performed utilising a series of hydrodynamic, cosmological N-Body simulations, studying the 
resulting power spectra and galaxy halo properties, such as density and temperature profiles. 

Results. In the cases of universal couplings, the deviations in the baryon fraction from ACDM are smaller than in the 
cases of non-universal couplings throughout the halos. The same is apparent in the power spectrum baryon bias, defined 
as the ratio of gas to dark matter power spectrum. Deviations of the density proHles and power spectra from the ACDM 
reference values can differ signihcantly between dark matter and gas because the dark matter deviations are mostly 
larger than the deviations in the gas. 


1. Introduction 

One of the most challenging problems in the field of cos¬ 
mology is to under stand the accelera ted expansion of the 
late-time Universe (Riess et al. 19981. The ACDM model 
is the most accepted explanation for this expansion and 
is reached by modifying general relativity (GR) through 
adding a cosmological constant, dark energy, to the energy- 
momentum tensor. 

An alternative to adding a dark energy component to 
GR is to modify gravity, by altering the Einstein-Hilbert 
Lagrangian, which the Einstein tensor is derived from. 


Many of these modified gravity theories exist (Brans & 
Dicke||1961[ |de Felice fc Tsujikawal |2010[ |Sotiriou| |2006t 


Clifton et al. 2012 Boehmer fc Mota |200 8 1 and several 
have previously been studied (|5rax et al. 2013 Barreira 


et al. 112013 Li et al. 2011[ U et al. 2012 Davis et al.||2012r 


gravity probes ([Bertotti et al. 11200*3 Will 


20141. 


If we assume that this fifth force acts on larger scales 
than the solar system, then some mechanism is needed to 
negate the fifth force on solar system scales. One way to 


found in the literature ( 

Khoury |2010 

Brax et al. 

2012 

Hinterbichler & Khouryl 

20101 IKhoury 

& Weltman 

2004 

Vainshtein 11972 

Koivisto et al.|2012l that screens the fifth 


20101, which screens the fifth force in regions of high den- 


Barrow fc Mota 2003[ Mota et a .| 2008 Puchwein et al.| 
2013 Winther et al. 2012[ |. These theories are often imple- 
mented by introducing a scalar field to the Einstein ten¬ 
sor, which couples to the matter component of the Uni¬ 
verse. This scalar field gives rise to a fifth force, an addi¬ 
tional gravitational force, which is negligible at solar sys- 
tem scales and below, according to labora tory experiments 
(]Hoyle et al. 2004| Dimopouloset_^ 2007 land solar system 


sity. 

With these modified gravity theories comes the chal- 
l enge of finding methods to test them ag ainst observations 
( Terukina et al.j 2014 IWilcox et al. 20151. Theorists in the 
past have mainly used predictions from models and simu¬ 
lations that only include dark matter due to the simplistic 
nature of dark matter, when constraining modified grav¬ 
ity theories. Use of models that only include dark matter 
can be justified by the matter composition of the Universe, 
which is 84.4 % dark matter and 15.6% baryonic matter, ac¬ 
cording to the Planck Gollaboration et al. (20151. However, 
astronomers observe the electromagnetic spectrum emitted 
from the baryonic components of the Universe, leaving the 
community with a disconnect between theories and observ¬ 
ables. 

Introducing the concept of a bias between the dark mat¬ 
ter and baryonic components is one way to rectify this dis¬ 
connect. The bias assumes that the behaviour of the two 
components are the same, but that the strength or am¬ 
plitude of the behaviour might be different. The commu¬ 
nity generally assumes that the bias of the standard model 
AGDM is equal to one, which justifies studying the Universe 
and comparing the observations to simulations that only 
include dark matter. However, if the bias is not unity re¬ 
searchers might greatly underestimate or overestimate their 
findings. 

The bias can deviate from unity even in AGDM because 
of baryonic effects other than a non-univers al coupling, such 
as AGN-feedback ( van Daalen et ar||2011 1 or the hydrody- 
namic effect observed in the Bullet cluster (|Tucker et al. 
1998 Glowe et al. 20061. This means that these baryonic 


we study the symmetron model (Hinterbichler & Khoury 


processes need to be comprehensively understood before the 
bias c an be used to test a non-universal coupling. Recent 
work (Schaller et al. 20151 has made progress in studying 
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the differences between galaxies (gas) and dark matter in 
the ACDM model. 

In our previous work ( Hammami et ar]|2015 ), we stud¬ 
ied the effects of adding a hydrodynamic gas to an exist- 
i ng N-body code wit h the symmetron model implemented 
(Llinares et al. 20141 with the assumption that the scalar 
field would couple to the gas and the dark matter to a 
universal coupling. There are, however, no justifications for 
assuming that the dark matter and gas have the same cou¬ 
pling to the scalar field. In fact, several known particles, 
such as neutrinos and baryons, do not interact in the same 
manner with electromagnetic forces, and it therefore stands 
to reason that it is worthwhile to study the effects of non- 
universal coupling between the scalar field and matter com¬ 
ponents. 

In this paper, we study density and temperature profiles 
and power spectra of dark matter and gas for symmetron 
models with different values for the coupling strength. The 
paper starts with an introduction to the symmetron model 
in Section]^ followed by a very brief section on the simula¬ 
tion parameters in Section [3| presenting and discussing the 
power spectra in Section ^and the density and tempera¬ 
ture profiles in Section and finishes with conclusions in 
Section |6] 

We work with some models that are extremely cou¬ 
pled to the scalar field to push any signatures of the non- 
universal coupling to their limits and possibly reveal sig¬ 
natures that would not immediately be clear from more 
sensible scalar field couplings. 


2. The symmetron model 


Introduced by Hinterbichler & Khoury (20101 the sym¬ 
metron model IS a scalar theory o f gr^ity using a sym¬ 
metric poten tial, where the action (Sotiriou 2006 Fujii & 
Maeda|[2003 1 is 


S= 


R ..2 1 




( 1 ) 


T fj,i> ^ ^ 


;1(Q R is the Ricci scalar, Mpi is the 
is the determinant of the metric 


where ijj is the scalar field 
Planck mass, and g = 
tensor in the Einstein frame, which can be converted to the 
Jordan frame by 




fiiy • 


( 2 ) 


The conformal factor satisfies A ~ 1 for the symmetron 
model and we use this approximation throughout. For more 
on these frames and the transformations between them and 


mami 


possible errors, see Faraoni et al. (19991 and Brown & Ham- 

(poT^. 


To preserve the behaviour of gravity, as described by 
GR, the symmetron model utilises a screening mechanism 
that triggers based on a set density value Passb a-t the solar 
system scale(a region of high density). In regions of low 
density, the symmetron would produce a modification of 


We only study the quasi-static limit (Llinares fc Mot^|2013 


Noller et al.|2014 \ of the scalar field, where time derivatives are 
ignored. 
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order one on the gravity. To accomplish this, the potential 
in the action above is defined to be symmetric, as in 

VW = Vo- ( 3 ) 

where is the scalar field, /i is a mass scale, and A a di¬ 
mensionless parameter. Likewise, the coupling factor is also 
symmetric, 

with M being another mass scale. 

To find the stress energy tensor for the symmetron 
model, we vary the action with respect to the metric 

r^. = A(^)TM+TW 

= ^('0) [{P + p)uit.u„ + Pg^iA (4) 


where P and p are the pressure and density, respectively. 

The scalar field component of the stress energy tensor 
is not covariantly conserved, i.e. 


^ 0 , 


while the total stress energy tensor is 


( Misner et al.|[l97^ 


= 0 . 


( 5 ) 


The equation of motion for the scalar field is found by 
varying the action again, this time with respect to the scalar 
field, 

□i/, = W(V') - (6) 


where is the trace of the stress energy tensor 

rp{m) ^ 

The right side of Eq. (|6 1 can be recognised as an effective 
potential, and, using Eq. (S, we write 


14ff(V’) = Vb + ^ + (7) 


With this potential, the scalar field goes to zero in re¬ 
gions of high density, pm S> while in regions of low 

density it reaches a minimum of The addition 

to gravity, the fifth force scales with the value of the scalar 
field, and we see from this that it is suppressed in regions 
of high density. 

We redefine the free parameters p,, M, and A to /3, Aq, 
and assB to a set of p arameters that are m ore physical 
intuitive as described in Winther et al. (20121, 


MpiV'o 

M2 


3 mlvtmMly 

«SSB - M2^2 


>2 — 
Aq — 


1 

V’ 


( 8 ) 

( 9 ) 

( 10 ) 


where Hq is the Hubble factor at present day (z=0) and 
rim is the matter density parameter. 
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The relative strength of the fifth force to the gravita¬ 
tional force is represented by /3, the moment of breaking 
symmetry is represented by the expansion factor assB = 

1 /3 

(flnioPco/psSB) ' , and the range of the fifth force is repre¬ 
sented by Aq in units of Mpc/h. 

A dimensionless scalar field y is defined as 



( 11 ) 


with an equation of motion in the quasi-static limit 
(Llinares et al. 


as 


^9 ^ 


/ aSSB Pm , 3 

- —X + X -X 

\ a J 


( 12 ) 


where is the mean density. 

The equation of motion for the position x of the dark 
matter N-Body particles has been derived from Eq. 0 in 
Llinares et al. (20141 and takes the form 


X + 2Hx + ^Vd) -I- 


1 A'U)„ , „ 


(13) 


Table 1. Coupling factor combinations explored. 


Configuration 

/3dm 

/^Gas 

DMIGI 

1.0 

1.0 

DMIOGIO 

10.0 

10.0 

DMO.lGO.l 

0.1 

0.1 

DMIOGI 

10.0 

1.0 

DMIGIO 

1.0 

10.0 

DM0.IGl 

0.1 

1.0 

DMlGO.l 

1.0 

0.1 


|e+v(up)=o, 

^ -H (D • V)u = -IvP - vl> - 

dt ^ ^ p Aii;) ^ 

8F P A'iilPi 

+ u • VA -b ^ • Vu = -(u-V)$-• VV'. 

at p A(V’) 


( 20 ) 

( 21 ) 

( 22 ) 


where $ is the Newtonian gravitational potential. 

The fluid equations for the symmetron model is a spe¬ 
cial case of the general fluid equations for a scalar-tensor 
theory. Via the action Eq. 0 . stress-energy tensor Eq. 
conservation law Eq. ([^, and working in the Newtonian 
Gauge, 

ds^ = -(1 -b 2<^)dP + 0^(1 - 2^)5ijdx^dx\ (14) 


the fluid equations are derived, 

dp 


— -b V(up) -bSiLp = 0, 

dv 1 

+ p) Hv + ^ + (t; • V)^; + 

I at . 

+vp + ^pVV' = o, 

dE P 

——b 2HE + V ■ VE H-• Vu 

dt p 


(15) 

(16) 

(17) 


where H = ^ is the Hubble factor, v is the velocity of the 
fluid, and E is the internal energy of the fluid. 

To remove explicit dependencies on a and H from the 
equations above, we use a variation of the super-comoving 
coordinates from Martel & Shapiro (19981, represented by 
a tilde. 


X = aX: di = a '^dt, p = Pp, v = a^v, (18) 

^ = a-ip, P = P'P, $ = 0^$, E = a^E] (19) 


all equations from this point on are in comoving coordi¬ 
nates. 

By excluding terms of second order and assuming static 
pressure, the field Eqs. (15 - 17) transform tc0 


^ Simulations beyond the static limit were presented in Llinares 
|& Mota| (|2013| |2014|), finding only sub-percent differences be- 
tween the static and non-static solutions. 

^ With this transformation, the derivative in A' {p) is now with 
respect to i/i. 


With this approach, the symmetron model version of 
the fifth force is 


f^ = -^M^vP = -- 

^ A{P) 


P 

M2 


P 




-2^P^-Jp^P 




2 


mAiQ 3 


xVx- 


(23) 


SSB 


For more on the symmetron model, see [Hinterbichler fc] 


Khoury (2010). 


3. Parameters 

The coupling factor defined above is split into two new cou¬ 
pling factors, 

, (24) 

tPGas 

which replace the coupling factor in the dark matter and 
fluid equations, respectively. 

In order to study the effect, we choose couplings of vary¬ 
ing orders of magnitud^ and our chosen configurations are 
shown in Table 1. 

The simulations were run using 1024 cores, 256^ dark 
matter particles, with a box width of 256 Mpc/h, and 
six levels of refinements. The background cosmology is a 
standard ACDM background, with h = 0.65, Ha = 0.65, 

The extreme couplings with /3 = 10 can induce accelerations 
upwards of F/j = 200 Fgr. Accelerations of this kind can result 
in relativistic velocities, requiring a relativistic set of equations 
to properly describe the systems. Luckily, none of our models 
induced relativistic velocities, the fastest dark matter particle 
in our simulations reached a speed of Umax = 0.033c for the 
DMIOGIO model, barely a factor of 2 higher than the fastest 
dark matter particle in the ACDM model with Umax = 0.016c. 
We continue with the extreme models to push any signatures of 
non-universal coupling to its limit. 
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flm = 0.35, and ft}, = 0.05. The chosen symmetron model 
has assB = 0.33 and Aq = 1 Mpc/h. 

Because of the use of extremely coupled models, we 
ran a set of convergence tests to verify that the errors in¬ 
duced by extreme coupling were not too extensive. The tests 
showed that the code handled the extreme models well for 
the most part, however, for the power spectra there were 
slight errors at the smallest scales fc > 4 Mpc/h and for the 
halo profiles at radius above R > 3i?2ooc- These errors were 
not significantly large, but results from the extreme models 
should be taken with a grain of salt in these regions. 

All results in the following sections only focus on the 
present day epoch, which corresponds to z = 0. 


4. Power spectra 


With the use of the open POWMES code (Colombi & 


Novikov 20111, we compute the power spectra for both dark 


matter and gas. To calculate the gas power spectrum, we 
treat each cell as a particle with a mass defined as 


m = pVceW, (25) 

where p and fyeii is the gas density and volume of the cell, 
respectively. 

In Fig.[T]we present the power spectra for all our models 
and the deviations of these power spectra from ACDM for 
both the dark matter and gas components. 

At the large scale range of the power spectra, k < 1 
Mpc/h, the extreme dark matter models (DMIOGIO and 
DMIOGI) show stronger effects of the scalar field in dark 
matter power spectra than in gas power spectra. This is 
because the power spectra of a component are most sensi¬ 
tive to changes to that component. The smaller differences 
between the various models in the gas spectra are due to 
the simulation, which is strongly dominated by dark mat¬ 
ter compared to the gas, so that the extreme gas model 
(DMIGIO) is suppressed by dark matter. This also accounts 
for that fact that DMIGIO is not clearly distinguishable in 
dark matter spectra. 

On smaller scales, k > 1 Mpc/h, the variations between 
the models are most noticeable in gas spectra, as is to be 
expected as the baryonic processes are strong at the smaller 
scales. Overall, we see a clear example that big deviations 
in dark matter power spectra do not mean big deviations 
in the gas power spectra. 


4.1. Power spectra deviations 

The main difference between DMIOGIO and DMIOGI in 
dark matter power spectra deviations is that the latter is 
shifted slightly towards smaller scales and has a slightly 
smaller amplitude. The addition of the extremely coupled 
gas increases clustering on larger scales, while decreasing 
clustering on smaller scales. 

The gas spectra deviations show that the spectra is 
greatly influenced by both the dark matter and gas. This 
is evident because DMIOGIO has the strongest deviations, 
while the second and third largest deviations come from 
DMIGIO and DMIOGI, respectively. The effect of dark 
matter on the gas spectra is much stronger than the ef¬ 
fect of the gas on the dark matter spectra, a consequence 
of the dominance of dark matter. 
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k (h/Mpc) 


Fig. 2. Deviation of the bias from the ACDM bias for our 
various models. The top image shows the extreme models, the 
bottom image the remaining models. All results are at z = 0. 


DMIGI and DMlGO.l show that a minimally coupled 
gas reduces the amplitude of the power spectra in the range 
k G [0.1,3), with the largest difference from DMlGl at 
fc ~ 1 Mpc/h. However, at the smallest scales these models 
end up with the same amplitude, indicating that at this 
point the gas coupling is irrelevant for dark matter spectra. 


The models DMO.lGO.l and DM0.IGl reveal that the 
gas has an impact on dark matter spectra, although not sig¬ 
nificantly, and the same very minor effect was also evident 
when comparing DMIGIO to DMIOGIO and DMIOGI. 


The gas spectra is more susceptible to changes to dark 
matter coupling, however, it is still most dependent on gas 
coupling, as demonstrated by the increasing amplitude as 
DM10G1<DM1G10<DM10G10. All the models with an ex¬ 
treme coupling have very low amplitude at small scales, and 
at the very smallest scales, they all end up at more or less 
the same value regardless of whether the coupling is univer¬ 
sal or not. The same effect is seen when comparing DMIGI 
to DMO.lGl. 


10 










































A. Hammami and D. F. Mota: Non-universal coupling in screened scalar-tensor gravity 





4 



3.5 


3 

S' 

2.5 

s 

Q 

o 

< 

CL 

2 

1.5 

S' 

1 

E 

CL 

0.5 

0 


-0.5 


-1 


k (h/Mpc) 



0.01 0.1 1 



k (h/Mpc) k (h/Mpc) 

Fig. 1. Top: power spectra for all our models. Middle: power spectra deviations of our extreme models from the ACDM model. 
Bottom: power spectra deviations for the remaining models. Left column shows the dark matter component, while the right column 
shows the gas component. All results are at z = 0. 


4.2. Bias 

The bias is defined as the ratio between the gas power spec¬ 
trum and the dark matter power spectrum, 


-^Gas 

fbM ’ 


(26) 


and is shown in Fig.[^ In the very large scale region, k < 0.3 
Mpc/h, the DMIGIO model has a bias that sky rockets 
to values larger than the ACDM bias. On the other hand, 
DMIOGI exhibits the exact opposite behaviour with a lower 
bias that rapidly decreases. DMIOGIO displays behaviour 
intermediate between these two models, but much closer to 
DMIOGI. 

The remaining large scales, 0.3 < fc < 1 Mpc/h, prove 
a turning point for all models. DMIGIO stops its growing 
bias and starts to decrease, and DMIOGI slows its rapid 
decrease and seems to stabilise at a constant value with 


= 0.25&ACDM- Similarly, DMIOGIO starts to level out, 
however, as DMIGIO starts to decrease DMIOGIO does as 
well. 

On smaller scales, the biases of these models start to 
diminish and eventually all are at lower values than the 
ACDM bias. Bias deviations in the smaller, non-linear 
scales are strongly correlated with gas coupling, while also 
dependent on the coupling to the dark matter. DMIOGIO 
reveals that the components compound the effect of the de¬ 
viations so that the DMIOGIO deviation is larger than the 
sum of the DMIOGI and DMIGIO deviations. DMO.lGO.l 
exhibits almost no deviations from ACDM, as is to be ex¬ 
pected. In the final models, the behaviour of the bias de¬ 
viations act in extremely different manners depending on 
what components are strongly coupled. 

For DMlGO.l, we see that the bias deviations plummet 
at large, non-linear scales until they reach more or less con¬ 
stant values that are decreasing slightly from /c ~ 1 Mpc/h 
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and out. The DM0.IGl model exlhbits the exact opposite 
behaviour at large, non-linear scales. DMIGI displays a be¬ 
haviour intermediate between the other two models at the 
large scales, and then eventually starts to decrease to much 
less than the sum of the respective deviations. The power 
spectrum bias at larger scales shows that models with a 
universal coupling have fewer deviations from the AGDM 
model than the models with a non-universal coupling. 

For observational astronomers the implications of this 
behaviour is that if the bias deviations are greater than 
unity, 6^ > 6acdm, then researchers who infer dark matter 
properties from baryonic physics features will make predic¬ 
tions with values of the power spectrum that are too high. 
The opposite is true if the bias deviations are less than 
unity, < 6acdm- 


5. Halo profiles 


In this section, we present density and temperature profiles 
for multiple halos identified by u sing the Rockstar code 
developed by Behroozi et al. (2013). We study massive halos 
with mass in the range [lx lO^'*/i“^M0,5x the 

high mass of these halos should ensure that the screening 
mechanism are triggered in dense regions. 

In previous works (Hammami et al.|2015 Llinares et al. 


20141, the halos that had not yet reached a relaxed state 
were filtered out following the me thods described in |Neto| 
et al. (20071 and Shaw et al. (20061. The method used rela¬ 
tions between the kinetic ana potential energy and surface 
pressure to determine if a halo was relaxed or not. |Gronke| 
et al. (20141 refined the method to take the effects of mod- 
itied gravity in the virialisation state of the halos into ac¬ 
count. However, if we filtered the non-relaxed halos, this 
would eliminate too many halos in the extreme models and 
leave almost no remaining halos. For this reason, we do not 
filter the non-relaxed halos from the relaxed halos. 

The halo profiles are calculated by sweeping over all 
the cells in the simulations, binning dark matter particles 
and baryonic gas properties in annular bins for each halaH 
then averaging over all halos. The profiles range from 10% 
of the virialisation radius, r = O.li? 200 cj to ten times the 
virialisation radius, r = 10i?2OOc- This range was chosen to 
properly catch all behaviours of the fifth force on the dark 
matter and gas halos, while also avoiding the inner regions 
of the halos where the resolution of our simulations is low. 
All profiles are calculated at the present epoch z = 0. 


5.1. Density profiles 

The density halo properties are presented in Fig. An 
extreme coupling in the dark matter strongly affects both 
dark matter and gas density profiles. For dark matter pro¬ 
files, the effect is mostly contained at the exterior of dark 
matter halos, while the gas profiles show signs of extreme 
dark matter coupling in all parts of the halo and outside. 
The diminished clustering at the inner regions of the ha¬ 
los can be explained with an environmental effect from the 
dark matter on the gas. 


® Rockstar filters out unbound particles when calculating the 
halo properties, while we construct the profiles by sweeping over 
all cells and binning particles within r = 10i?200c, potentially 
including these unbound elements. 
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The dark matter clusters faster than the gas because of 
its collisionless nature, in which the gas is prevented from 
collapsing, due to the pressure, and the dark matter clus¬ 
ters unhindered. This means that the dark matter reaches 
higher densities at a faster rate than the gas, and from the 
description of the screening mechanism, we know that the 
screening is triggered by a combined density, dark matter 
plus gas, threshold, resulting in the dark matter triggering 
the screening mechanism before the scalar field has had a 
chance to work on the gas component as much as on the 
dark matter. 

An extremely coupled gas has a minor, however, not 
negligible effect on the dark matter profiles, while it has a 
huge impact on the gas profiles. This is expected as there is 
much more dark matter than gas, resulting in the effect of 
the gas not being as strong. Also, DMIOGIO is the model 
with the biggest effect on the gas profiles, while the dark 
matter profiles seem to prefer DMIOGI. 


5.2. Density profiles deviations 

The deviations from AGDM for the halos are found in the 
lower two rows of Fig. The deviations confirm the conclu¬ 
sions from the total density profiles, and also reveal previ¬ 
ously unseen effects. All models without extremely coupled 
dark matter make dark matter halos cluster less at the inner 
regions than in those two cases (DMIOGIO and DMIOGI). 
The exact opposite occurs in the gas power spectra, where 
all extremely coupled models (the two previously mentioned 
and DMIGIO) cluster less than the other models. This is 
due to baryonic physics preventing the gas from collapsing 
as far inwards as the dark matter, and then the gravita¬ 
tional contribution from the gas on the dark matter also 
prevents the dark matter from clustering. 

The peak of deviations does not occur at the same radius 
for all the models. For DMIGIO, the peak in deviations is 
closer to the halo centre, while for DMIOGIO and DMIOGI 
the deviations peak further out. This effect is due to the ex¬ 
treme coupling in the gas. The extreme coupling allows the 
scalar field to affect the gas component to such an extent 
that it collapses further inwards before the environmental 
effect of the combined density triggers the screening mech¬ 
anism. 

The bottom row reveals that DMIGIO has a strong ef¬ 
fect on dark matter halos throughout the density profiles 
with higher peak deviation than DMlGl. The gas profiles 
are under-dense at the inner regions of the halos in the cases 
where the dark matter is normally (or extremely) coupled to 
the scalar field. In the cases where the dark matter is min¬ 
imally coupled, the profiles approach AGDM at the centre 
of the halos. This behaviour further asserts our conclusions 
about the environmental effect. 

To study how deviations from AGDM differ in dark mat¬ 
ter and gas cases, we introduce the deviation bias (5^^, de¬ 
fined as the relative difference between the deviations 


jDM _ ^DM ~ Acas 
^Gas 


Ppm —PACDM _ PGae —PAGas 

PACDM_PAGas 

PGaa —PAGaa 

PAGaa 


(27) 


where pAGas is the gas AGDM density and pacdm is the 
DM AGDM density. 

The DMIOGI model represents a model that has an 
enormous deviation from AGDM in the dark matter com¬ 
ponent of Adm ~ 8.5 at i? ^ 4i?200c) and a much smaller 
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Fig. 3. Density properties for the massive halos with mass in the range [lx x Top: density profiles for 

ACDM and symmetron models. Middle: density profile deviations from ACDM for the extreme models. Bottom: density profile 
deviations from ACDM for the remaining models. Left column shows the dark matter component, while the right column shows 
the gas component. All results are at z = 0. 


deviation of Acas 2.5 in the gas component, giving a de¬ 
viation bias of « 2.4. The same effect, but now with 
enourmous deviations in the gas component, is present in 
DMIGIO where the deviations from AGDM are Acas 1-5 
for the gas and Adm ~ 0.5 for the dark matter, giving a 
deviation ratio of ~ — I ■ 


5.2.1. Baryon fraction profiles 


One of our aims is to find a method of distinguishing be¬ 
tween models that have universal coupling and models with 
non-uni v ersal coupling. For thi s purpose, we use the baryon 
fraction ( Kravtsov et al.|[2005 |, 

/Gas = - ^ - , (28) 


PGa 


Pdm 


and present the baryon fraction profiles in Fig. 

The baryon fraction shows that all the models have 
the same behaviour as AGDM, except the three models 


DMIOGIO, DMIGIO, and DMIOGI. DMIOGIO converges 
with the AGDM baryon fraction earlier than DMIGIO and 
DMIOGI, however, in general this model has a higher devia¬ 
tion at the inner region. In fact, the models with a universal 
coupling seem to deviate less from AGDM than the models 
with a non-universal coupling. 

This is the same behaviour displayed by the power 
spectrum bias at larger scales, and allows us to conclude 
that significant deviations in the baryon fraction from the 
AGDM model throughout the halos might be an indication 
of a non-universal coupling. 


5.3. Temperature 


The temperature is a very interesting component to study 


due to its close rela tion observables (Wilcox et al. 2015 
Terukina et al.|[^014l. The temperature is not an output of 


our code and needs to be reconstructed using the ideal gas 
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Fig. 4. Baryon fraction profiles for the halos with mass in the 
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p = RspT, (29) 
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where p is the thermal pressure, i?s = ^ is the specific 
gas constant, m = 0.59m// is the mean mass of the gas, 
ttih is the hydrogen mass, and p is the gas density. The 
temperature profiles are made and analysed in the exact 
same manner as the density profiles. 


Fig. 5. Temperature properties for the our halos. Top: tem¬ 
perature profiles for ACDM and symmetron models. Middle: 
temperature profile deviations from ACDM. Bottom: deviations 
with the extremely coupled models filtered out. All results are 
at 2 = 0. 


5.3.1. Total temperature profiles 

The temperature profiles are presented in Fig. The tem¬ 
perature is a product of the baryonic processes and is most 
sensitive to strong coup ling between the gas and scalar field 
( Hammami et al)]|2015 |. 

The stronger the coupling to the scalar field is, the 
higher the temperature in the halos is. The slope of the 
profiles outside of the halo, as the temperature starts drop¬ 
ping, is shallower for DMIOGIO, DMIGIO, and DMIOGI 
than in the other models. DMIOGIO gives an even higher 
temperature than DMIGIO, indicating that the dark mat¬ 
ter plays a vital role in the temperature of baryonic halos. 
This is also evident from DMIOGI, which has a higher tem¬ 
perature than the less extreme cases. 
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5.3.2. Temperature profiles deviations 

The temperature deviations from the AGDM temperature 
are found in the bottom two rows of Fig. DM0.IGl and 
DMO.lGO.l are two similar models in which one has uni¬ 
versal coupling and the other does not. These two models 
display no signatures that can be interpreted as a trait of 
non-universal coupling. DMlGO.l and DMO.lGO.l, on the 
other hand, display an increase in temperature outside of 
the virialisation radius. 

Both of these effects come from the fact that these 
halos are so massive that they are most likely screened 
from the scalar field. The difference between DMlGO.l and 
DMO.lGO.l is that the regions outside of the halos are non- 
screened, and the stronger coupling to dark matter allows 
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DMlGO.l to be influenced by the scalar fleld to a larger 
extent than DMO.lGO.l. 


6. Conclusions 

We investigate scalar-tensor theories of gravity, which 
present a non-universal coupling. That is, baryons and 
dark matter are coupled with the scalar degree of freedom 
with different coupling strengths. The models we investigate 
utilise a screening mechanism to suppress the deviations 
from GR at small (solar system) and large cosmological 
scales: the symmetron screening mechanism, specifically. 

As a result of the screening mechanism, the strongest 
signatures in these models are expected to occur at the non¬ 
linear regime of structure formation. Therefore, in order to 
unveil the imprints of these theories at astrophysical scales, 
we ran several hydrodynamic cosmological N-body simula¬ 
tions. We compared models with and without a universal 
coupling to the symmetron scalar fleld, and showed that 
several astrophysical observables (density profiles, temper¬ 
ature profiles and power spectra) show significant differ¬ 
ences between the dark matter and gas components when 
the coupling is non-universal. 

The deviations from AGDM are typically larger in the 
gas than in the dark matter near the centre of the halo. The 
opposite holds true at larger radii, where the dark matter 
deviates more strongly from AGDM. However, this is not 
the case for models in which coupling of the gas is signifi¬ 
cantly stronger than that of the dark matter. 

For power spectra, dark matter deviations are larger 
than that of the gas in models with universal coupling or 
in models in which dark matter coupling is stronger than 
gas coupling. 

Our attempt to And signatures in density profiles and 
power spectra, which would reveal whether coupling to the 
scalar fleld is universal or not, revealed one signature: in 
the cases of universal coupling, the deviations in the baryon 
fraction and bias from AGDM are smaller than in the cases 
of non-universal couplings throughout the halos. This is 
expected, since GR is a universally coupled theory. If ob¬ 
servers And the baryon fraction or power spectrum bias to 
deviate from the calculated AGDM bias, then this might 
very well be a sign of a non-universal coupling, and there¬ 
fore a breaking of the equivalence principle. 

Separating the dark matter and gas will prove to be a 
challenge for observers intending to compare their results 
with theories, since the dark matter is not a direct observ¬ 
able, while the gas is observable. To work around this, we 
propose to use the baryon fraction as presented in this pa¬ 
per. To measure this baryon fraction, first we suggest the 
constructing of the total density profile, pxot = Pdm + Pgas) 
using rotational velocity profiles of individual galaxies to 
measure the dynamical mass, and then constructing the 
gas density profile Pgasi using X-ray surface brightness and 
X-ray temperature profiles. 

A caveat with the above mentioned method is that the 
calculated dynamical mass is dependent on the particular 
gravity model used. An alternative option is to calculate the 
lensing mass using gravitational lenses as the path of light 
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The above method for detecting non-universal couplings 
may not be possible with the current state of observa¬ 
tional and theoretical limits. The deviations we found from 
AGDM are all quite small, when we exclude the extreme 
models, and also contain large uncertainties depending on 
the modelling of feedback physics. 

Furthermore the observational baryon census of galaxies 
contains significant errors and ma y still be incomple te, i.e. 
the halo missing baryon problem ( Werk et al.| 20141, while 
the baryon consensus in galaxy clusters is less challenging as 
the majority of the gas is hot enough to be visible in X-rays. 
The above method is therefore more likely, yet still very 
challenging, to work in galaxy structures than in galaxies. 

If one takes seriously the possibility of matter compo¬ 
nents with a non-universal coupling to a gravity scalar de¬ 
gree of freedom, then our work shows the bias will be greatly 
affected. Therefore, attempts to rule out or constrain mod¬ 
ified gravity theories by comparing dark matter predictions 
to the observed quantities based on baryonic properties may 
be misleading, and one must consider the possibility of a 
non-universal coupling that might skew the conclusions and 
dark matter properties that are inferred from baryons. 
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